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We theoretically study the properties of Nambu-Goldstone bosons in an interacting single- 
component Bose-Einstein condensate (BEC). We first point out that the proofs of Goldstone's 
theorem by Goldstone et al. [Phys. Rev. 127 (1962) 965] may be relevant to distinct massless modes 
of the BEC: whereas the first proof deals with the poles of the single-particle Green's function G, 
the second one concerns those of the two-particle Green's function. Thus, there may be multiple 
Nambu-Goldstone bosons even in the single-component BEC with broken (7(1) symmetry. The 
second mode turns out to have an infinite lifetime in the long-wavelength limit in agreement with 
the conventional viewpoint. In contrast, the first mode from G, i.e., the Bogoliubov mode in the 
weak-coupling regime, is shown to be a "bubbling" mode fluctuating temporally out of and back into 
the condensate. The substantial lifetime originates from an "improper" structure of the self-energy 
inherent in the BEC, which has been overlooked so far and will be elucidated here, and removes 
various infrared divergences pointed out previously. 



I. INTRODUCTION 

On the basis of the exact equations of motion for the 
single-particle Green's function G and condensate wave 
function, 1 , it was shown in a previous paper— that the 
poles of the two-particle Green's function JC in an in- 
teracting Bose-Einstein condensate (BEC) are generally 
not shared with those of G. This implies that there 
may be multiple excitation branches even in a single- 
component BEC, contrary to the widely accepted view- 
point of attributing a unique phonon branch due to den- 
sity fluctuations.— ~— Hence, it is desirable to clarify the 
basic properties of the poles of G and JC. Here, we will 
perform it in terms of Goldstone's theoremSr— based 
on a self-consistent perturbation expansion that simul- 
taneously satisfies Goldstone's theorem and conservation 
laws. 1 

Goldstone's theorem on spontaneously broken contin- 
uous symmetry has crucial importance over a wide range 
of fields in modern theoretical physics. According to a re- 
cent review fii it is considered to predict a single massless 
mode called the Nambu-Goldstone boson for the single- 
component BEC with broken U(l) symmetry. Since the 
mode is expected to dominate the thermodynamic prop- 
erties of the BEC, clarifying the basic properties of the 
Nambu-Goldstone boson(s) forms an essential prerequi- 
site to any theoretical investigations of the BEC. In this 
connection, it should be pointed out that basic proofs 
of Goldstone's theorem were given by Goldstone, Salam, 
and Weinberg (GSW) mainly in two different ways^ We 
will show that they may be relevant to distinct massless 
modes in general, contrary to the standard treatmen t 10 ^ 1 
where no explicit distinction between the proofs seems to 
have been made. 

The main results of the present paper are summarized 
as follows. For the single-component BEC, the two proofs 
by GSW ? are relevant to the poles of different Green's 
functions. The first one deals with the single-particle 



Green's function and can be identified with the proof 
of the Hugenholtz-Pines theorem^ The corresponding 
massless mode reduces to the Bogoliubov mode in the 
weak-coupling regime^ which has been regarded as a 
propagating mode with an infinite lifetime in the long- 
wavelength limit based on the work by Beliaevjii see also 
refs. [HI and [l6| on this point. We will show, however, 
that it is essentially a "bubbling" mode with a short 
lifetime, i.e., a mode that fluctuates temporally out of 
and back into the condensate, originating from an "im- 
proper" structure of the self-energy inherent in the BEC. 
To put it another way, the first mode of the BEC is not 
a well-defined "quasi-particle" ( "particle" ) in the termi- 
nology of condensed matter (particle) physics. Both the 
improper structure of the self-energy and the resultant 
bubbling feature of the first mode have escaped consider- 
ation so far and will be studied in detail. They also have 
the effect of removing infrared singularities pointed out 
by Gavoret and Nozieresj^ hence JC, for example, can be 
calculated safely without any artificial procedures. The 
second proof concerns a correlation function of three field 
operators. Its poles will be shown (i) to be common with 
those of JC and (ii) to have a vanishing imaginary part 
in the long-wavelength limit. In other words, the second 
mode is an observable propagating one. Thus, the two 
modes in the BEC from the two proofs are completely 
different in character. 



We now derive and elaborate on the above results. In 
§2, we recapitulate the main results of refs. G]and[2 

that 

are relevant to the present study. In §3, we consider 
the implications of Goldstone's theorem for the single- 
component BEC. Section 4 forms the core of the present 
study, in which we clarify the basic properties of the poles 
of G and JC. Section 5 provides concluding remarks. We 
set h = &B = 1 throughout with /cb the Boltzmann con- 
stant. 
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II. SUMMARY OF PREVIOUS RESULTS 

The system we consider consists of identical particles 
with mass m and spin described by the Hamiltonian 



with 



H — Hq + H ln t , 



H = J dVx^rOi^ri), 



(1) 



(2a) 



d'V / dV^^t^y^-rXr'j^r). 



(2b) 

Here ip^ and ip are the creation and annihilation opera- 
tors of the Bose field, respectively, K\ = — fi 2 W\/2m — fi 
with (j, the chemical potential, and V is the interaction 
potential. Though disregarded here, the effect of a trap 
potential may be included easily in K\ . Let us introduce 
the Heisenberg representations of the field operators in 
the Matsubara formalism as 17 




,TlH 



V> f (ri)e 



-TlH 

-nH 



(3) 



where the argument 1 in the round brackets is defined by 
1 = (ri,n), and the variable t\ lies in < t\ < f3 = 1/T 
with T denoting the temperature. The operators V'i(l) 
and ^2(1) were written as and ^(l) in ref. [l], respec- 
tively. Distinguishing them by the subscript i — 1,2 as 
above enables us to simplify the notation considerably, 
as seen below. 

The operator ?/>i(l) acquires a finite grand-canonical 
average in the condensed phase. Hence, it may be written 
as 

^ i (l) = * i (l) + i (l), (4) 
where ^(l) is the condensate wave function 

rry p -/3fl" ,/ ..(1 ^ 

ee ir ^ e _g ij ee (Ml)) (5a) 

with Tr the trace, and (<pi(l)) = by definition. Note 
that *i(l) = ^(1) = *( r i) m equilibrium with * de- 
noting the complex conjugate. Using 0^(1), we then in- 
troduce Green's functions as 

G y (l,2) = -(T T ^(l)03-,(2))(-l)'''- 1 , (5b) 

with T T the "time" -ordering operator.— They satisfy^ 

G u (l,2) = (-ir+^ 1 G 3 _ J - 3 - l (2,l) 

= (-l) l+ ''G*,(r 2 T 1 ,r 1 T 2 ). (6) 

The elements Gn and G12 were denoted as G and 
F in ref. Q], respectively, following the convention of 
superconductivity*^ 

We now summarize the main results of refs. [l] and 
that are relevant to the present consideration. 



A. Hugenholtz-Pines relation and Dyson-Beliaev 
equation 

Let us introduce the 2x1 vector ^ = [*i ^ 2 ] T and 2x2 
Nambu matrix G = (Gy), where T denotes transpose. 
They obey^ the generalized Gross-Pitaevskii equation (or 
Hugenholtz-Pines relation) 

G- 1 (l,2)a 3 *(2)=0 (7a) 

and the Dyson-Beliaev equatio n 17 i 18 

G- 1 (l,3)G(3,2)=a 5(l,2), (7b) 

respectively, where integrations over barred arguments 
are implied. The quantity G _1 is defined by 



G '(1.2)= (-a -£--a 3 K 1 



*(l,2)-f)(l,2), (8) 



where &o an d CT3 are the 2x2 unit and third Pauli ma- 
trices, respectively, K\ is given after eq. ([2]), (5(1,2) = 
S(ti — r 2 )<5(r 1 — r 2 ), and E denotes the self-energy ma- 
trix. 

Equation (ITal) was derived in Appendix A. 4 of ref. 
1 by a formal argument using the global gauge trans- 
formation relevant to the broken U(l) symmetry. In 
particular, the proof is free from a couple of assump- 
tions by Hugenholtz and Pinesi2 that (i) the system 
is homogeneous and (i) the self-energy is proper. The 
first row of eq. (|7a|) in equilibrium is written explic- 
itly as -Ki*(rx) = E u (l,2)tf(f 2 ) - E 12 (l, 2)**(r 2 ). 
By approximating En (1,2) = 2gS(l, 2)|*(ri)| 2 and 
E 12 (l, 2) = g5(l, 2)[1'(r 1 )] 2 for V(n - r 2 ) = gSfa ~ r 2 ), 
it reduces to the standard Gross-Pitaevskii equation^ - — 



Setting W — > y/n and K 



-fi in the same equa- 



tion with no the condensate density, we also obtain the 
Hugenholtz-Pines relation for homogeneous systemsji^ 
It was showni that the elements of E are obtained from 
a single functional $ by 

(5$ 



E y (l,2) = -2/3 



*Gji(2,l)' 



(9a) 



where $ = $[G, F, F, * 2 ] with G(l,2) = [G u (l,2)- 
G 22 (2,l)]/2, F(l,2) = [G 12 (l,2) + G 12 (2,l)]/2, and 
F(l,2) = -[G 2i (l,2) + G 2 i(2,l)]/2. Here, we have in- 
corporated the first symmetry of eq. manifestly in 
$; this procedure is indispensable for discussing the two- 
particle irreducible vertices in terms of It follows from 
eq. (|9a|) that Ey also obey eq. (j6]) for dj. There is an- 
other identity in terms of which is given by eq. (21b) 
of ref. 0], i-e., 

P 5*7^(1) = ^(MX-l)?- 1 *?®. (9b) 

This is the relation that makes eq. (|7al) a station- 
arity condition of the thermodynamic potential 57 = 
— lnTre - ^^ with respect to a variation of 

For a given eqs. © and (|5aj) form a closed set of 
equations that can be used to determine \& and G self- 
consistently. 
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B. Interaction energy 

By comparing eq. (O with the equations of motion for 
tfri (1), it was shown that the thermodynamic average of 
eq. (j2bf can be expressed rigorously in terms of ^, G, 
and E as eq. (12) of Ref. 1. In the present notation, it 
reads 

(H int ) = isjjCl^JK-iy-^^^lJ-G^S,!)]. 

(10) 

In §IV A| we will clarify a peculiar structure of the self- 
energy in the BEC based on this expression. 



C. Higher-order Green's functions 

We now consider a couple of higher-order Green's func- 
tions defined by 

(li|£|32jy) = -(T r V^(l)V' J (2)V'3-fe(3))(-l) fe - 1 

+* 4 (l)(T T ^(2)^ 3 - fc (3))(-l) fc - 1 , (11a) 

(12y|/C|43, fc ) = (T T V>i(l)^_j(2)^ fc (3)^-i(4))(-l) 3 ' +J 
-<T T ^(l)^3-,-(2))<r T V»fc(3)^3_i(4)) 
x(-l) i+ '. (lib) 

The two-particle Green's function JC can describe various 
collective modes such as density fluctuations. To derive 
equations of motion for £ and /C, let us add an extra 
perturbation given by the S matrix 



eq. (jllbl) is connected with eq. (113b[) as (12,-j|/C|43;fc/ 



S(P) = T T exp 



1 



^(i^-^x-ir 1 ^,!) 



The condensate wave function and single-particle Green's 
function in the presence of the nonlocal potential U = 
(Uij) are defined byii 



^(1,2) 



t, (T T 5(/3)^(1)) 
1 1 J " (508)) ' 

(T T 5W i (l)V. 3 _ J -(2))(-l)J- 1 

W)) 

Gg(l,2)-*P(l)*g_ i (2)(-l)^- 1 , 



(13a) 



(13b) 



where we have used the correspondent of eq. (j4|) in the 
second equality of eq. (|13b[) . Quantities ^(1) and 
G^(l,2) reduce to eqs. flU} and (|5bJ as U -> 6, re- 
spectively. It may be seen easily that £ of eq. (|lla[l is 
obtained from eq. (|13a[) by 



(1<|£|32 W ) = 2 



'^•(3,2)' 



(14a) 



where the limit C/ — 5- is implied after the differen- 
tiation; we will use this convention below. Similarly, 



2J0#(l,2)/fl7 Ifc (4,3) J Le. 



5GE(1,2) 



[* i (l)(2 3 _ J -|£|43, fc ) 



+ (l l |£|43 ife )* 3 -,(2)](-l)^ 1 . (14b) 

Thus, we only need to consider linear variations of ^E^ 7 
and Gfj to obtain C and K explicitly. 

To carry this out, we start from eq. (JTJ). The pertur- 
bation given by eq. (|12p produces in the right-hand side 
of eq. © an extra term -U'(l,2) with^^^ 

E , (lt2)s ^(l,2 ) + (-l ) ^,_ A ,_ < (2,l )< 

2 

Varying J7 — > £/ + 5U and subsequently setting U = in 
the resultant eq. (O, we arrive at the first-order equations 



G-^l, 2)& 3 6$ u (2) = [6U'(1, 2) + <5E U (1, 2)] <j 3 *(2), 

(16a) 

(^(M^G^M) = [,W'(1,3) + 6E U (1,3)]G(3,2). 

(16b) 

We next express the first-order self-energy 8YF in terms 
of 8G U and S^ u as 

^.(l,2) = -2(12 u |r( 4 )|43 rfe )<5G^(4,3) 

+2(12 y |r( 3 )|%)(-l)^ 1 ^(3), (17) 

with 

, m , v l<5Ew(l,2) <S 2 $ 



2 5G /fc (4,3) "^(2,1)^(4,3)' 

(18a) 

= 2(-l)"'+ r {12, : ,|r< 4 >|43 u _ t .)*,-(3), (18b) 

where eq. ([9]) has been used to derive the second expres- 
sions. We also introduce a set of matrices 

(l 4 |f (3) |32,,) = 2(-l)'- 1 (l%|r( 4 )|32 fej )* r (4) 

= (-iy(23 jk \T^\l 3 ^), (18c) 
(h\t^\2 3 ) = 2(-l) S - 1 (13 l ,|r( 3 )|2 J )%(3), (18d) 
(12 y |x (0) |43 ife ) = G ii (l,4)G fcj (3,2) + (-l) fe +'- 1 

xG v3 - fc (l,3)G3-;,,(4,2), (18e) 
(12 l3 |^ (3) |3 fc ) = (-1)^^(1)4,3-^(3,2) 

+a w J(3,l)*3-j(2)], (18f) 
(l i |£ (3) |32 fci ) = (-l)-'"- 1 [^(5(1,3)^(2) 

+<y il 3-^(i,2)* 3 _ fc (3)] 1 (i8 g ) 



and vectors 



(12^1^ = ^(1,2), (1^1^ = ^(1,2). (19) 



The quantities in eqs. ()18a|) - f)18dj) form "irreducible" 
vertices of our BEC, as seen below in eq. ([20)) . 
Let us substitute eqs. (H"5l) and ([T7| into eq. (fT6) 
and make use of eqs. (I18cp and (|18d[) as well as 
the symmetry (12y|r (4) |43 /fc ) = (U kl \T^\2l, 



(-l) i +J- 1 (2l3- J ,3- l |r (4, |43 ifc ) originating from eqs. © 
and (I18a|) . We thereby obtain a closed set of equations 
for 5$ u and SG U as£ 



^G± i3) SU + GtWa 3 S$ u - Gt (3) SG u : 



SG U = lx (Q) SU + X (0) rWa 3 S$ u - x {0) L^SG u . 



2- 



By considering eq. (|14p . they can be expressed equiva- 
lently in terms of C and K = 28G U / SU as 



er 3 £ 




"g* (3) ~ 




G 6 ~ 




~f(2) 


_ £ (3)- 




a 3 t 


X . 






+ 


x (0) 




p(3) 


_r (4) 







(20) 

which may be regarded as the Bethe-Salpeter equation^ 
for the condensed Bose system and can be solved formally 
with respect to t and K/ £ Using eqs. (|14bl) and (|18fl) . 
we finally arrive at the explicit expressions^ 



£^-3x (c) (i (3) -£ (3 V 4) )^ 

K=yW + y^I^y^ 3 ' + * (3)*(c)f ( V 4) 



(21a) 



_£(3)^(c)^ 3 ) ; 



(21b) 



where x^ 4 \ X^ c \ ancl X are defined by 

xW = ^)-i + rW)- 1 , (22a) 
x (c) = - f( 2 ) + f (3) xWr( 3 )y 1 , (22b) 

x (4)-i +j: (3) (d- 1 -f^Y't^ 1 



x. 



(q) = 



(22c) 



with superscripts c and q denoting "condensate" and 
"quasi-particle," respectively. It is worth pointing out 
that the poles of x^ m eq. (|2"Tj) are cancelled by those 
of x m the denominator of eq. (|22b[) . 

We observe from eq. (|2"Tj) that the function x^ c \ which 
is characteristic of the BEC, appears in both £ and JC. 
This x^ has been identified as G in the diagrammatic 
analyses of K for homogeneous systems.— ~— However, x^ c ^ 
of eq. (|22bl) is clearly different from G due to the addi- 
tional contribution — f ^ + 1^^ x^H^ in the denomina- 
tor. The properties of the poles of G and x > especially 
their differences, will be one of the main topics below. 




00 



(a) d>" 




FIG. 1: Diagrammatic representations of (a) "I 1 ' 1 ' and (b) 
$' 2 '. A filled circle denotes r' ', and a line with an arrow 
(two arrows) represents G (either F or F) in eq. (I2.':>p as in 
the theory of superconductivity.— The number below each 
diagram indicates its relative weight. 



D. The functional $ 

The key quantity in the present formalism is the func- 
tional $ = 4>[G, F, F, \J/ 2 ] with 

G(l, 2) = [Gn(l, 2) - G 22 (2, l)]/2, (23a) 
F(l, 2) = [G 12 (l, 2) + G 12 (2, l)]/2, (23b) 
F(l,2) = -[G 21 (l,2) + G 21 (2,l)]/2. (23c) 

Indeed, the self-energies and irreducible vertices are de- 
termined from this $ by eqs. (Pal) and (|18a[) - (|18dp . re- 
spectively; using them, we can obtain (\E ( , G) and (£, JC) 
self-consistently by eqs. (O and (|20|) . respectively. More 
generally, choosing an approximate <J> amounts to fix- 
ing the whole Bogoliubov-Born-Green-Kirkwood-Yvon 
(BBGKY) hierarchy's 

It was shown 1 - that <& can be constructed order by order 
perturbatively in terms of the symmetrized bare vertex 

T^{1\',22') = V{v 1 -v 2 )5{t 1 -t 2 ) 

x[6(l,l')S(2,2')+6(l,2')6(2,l% (24) 

so as to satisfy eqs. ((9]) and (fT0|) simultaneously. The 
first- and second-order contributions are given diagram- 
matically in Figs. [Ha) and[ljb), respectively. The corre- 
sponding analytic expressions read 

$W _ Ar( )(ll',22')[2G(I,r)G(2,2') 

-F(l, 2)F(l\ 2') - 4G(1, l')*i(2)* 2 (2') 
+F(1, 2> 2 (i')* 2 (2') + F(l', 2')*i(l)*i(2) 
+* 2 (I')*2(2')*i(2)*i(l)], (25a) 
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$(2) = __Lr( )(li',22')r( )(33',44') 

°P 

x {G(2, 3')G(3, 2')G(T, 4')G(4, l') 

-2F(2', 3')F(2, 3)G(1,4')G(4, 1') 

+F(2\ 3')F(2, 3)F(I', 4')F(I, 4) 

-4G(2, 3')G(3, 2')G(I, 4')*i(l)* 2 (I') 

+4F(2', 3')F(2, 3)G(l,2')*i(4)*a(l') 

+2[F(2', 3')*i(2)*i(3) + F(2, 3)# 2 (2')* 2 (3')] 

xG(l,4')G(4,l') 

-2[F(2', 3')*i(2)*i(3) + F(2, 3)tf 2 (2')* 2 (3')] 
xF(l',4')F(I,4) 

+2G(2, 3')G(3, 2 , )*i(I)* 2 (4 / )*i(4)* 2 (r) 

+2G(2, 3')*i(3)* 2 (2')G(l, 4')*i(4> 2 (l') 

-4[F(2', 3')*i(2)*i(3) + F(2, 3)* 2 (2')* 2 (3')] 

xG(l,4')*i(4)* 2 (l') 

+F(2, 3)* a (3')* 2(2')F(l, 2)* a (3')*a(I') 

+F(2',3')*i(3)*i(2)F(l',4')*i(4)* 1 (l) 

+2F(2', 3')F(2, 3)*i(l)* 2 (4')*i(4)*a(l')}- 

(25b) 



III. GOLDSTONE'S THEOREM FOR 
SINGLE-COMPONENT BEC 

We now discuss Goldstone's theorem^— in terms of 
the single-component BEC. Bose-Einstein condensation 
is a textbook example of broken U(l) symmetry de- 
scribed by the Mexican hat potential 2 ^ as well as a proto- 
type of the off-diagonal long-range order i 28 i 29 These facts 
may make the single-component BEC one of the most 
desirable grounds for studying the implications of Gold- 
stone's theorem. It will be shown that the basic proofs of 
Goldstone's theorem by GSW— are relevant to the poles 
of distinct Green's functions for the BEC. 

The first proof concerns the single-particle Green's 
function G defined by eq. (|5b[) , which obeys the Dyson- 
Beliaev equation (|7b[k our G^ corresponds to the com- 
plete propagator A£ • of GSW* 9 - Now, the conclusion of 
the first proof can be stated in the present context as 
follows: the operator G" 1 in eq. (|7b[) should also deter- 
mine the condensate wave function W — [Wi V E' 2 ] T by eq. 
(fTaj) . Equation (fTaj) . i.e., the Hugenholtz-Pines relation 
for general inhomogeneous systems ?i may be identified 
with eq. (18) of GSW, which manifestly tells us that poles 
of G contain a "massless" mode. It should be noted, how- 
ever, that eq. (fTJ does not say anything about the nature 
of the mode, which will be elucidated in ^IVl 

We next focus on the second proof by GSW involving 
the use of a commutation relation. 9 Although it is appli- 
cable only to homogeneous systems, it has been referred 
to more frequently as the standard proof of Goldstone's 
theoremj 10 ! 11 In the context of the BEC, the proof is rel- 



evant to 

Pn dP^) = f d 3 r f°° dte^-"^<[n(r,i),^(0)]>, 

J J — OO 

(26) 

where n(r, t) = e lHt ip' (r)-0(r)e tHt is the density opera- 
tor in the Heisenberg representation, and we have chosen 
the argument of ip arbitrarily at the coordinate origin 
based on the translational symmetry. Setting p = in 
eq. yields 

/OO 
dte iut {[e im Ne- im ,^(0)]), 
-OO 

where N = J d 3 r t/)' (r)ip(r) is the particle-number op- 
erator, i.e., the generator of the gauge transformation 
relevant to the broken U(l) symmetry. Since [N,H] = 
and [N,ip(0)] = —ip(0), we immediately obtain 

p„</,(0, oj) = -2TTy/n^S(uj), (27) 

which is finite for y/no = (rpi) ^ with a sharp peak 
at w = 0. Assuming that p n ^{p,uS) is continuous for 
(p,w) — > (0,0), which is justified here) 10 ' 11 we conclude 
that p„^,(p,a;) contains a massless mode with an infinite 
lifetime in the long- wavelength limit. 

It is clear from eq. (f2"6"j) that the second proof deals 
with a correlation function of three field operators, whose 
poles may be different from those of eq. (|5b[) . This is 
confirmed by relating p n ^ to C of eq. (|lla[) . Consider 
the Fourier transform 

£ni(p,^)= /d 3 r 21 y%r 21 (l 1 |£|22 11 )e^ p - r21+z ^ 21 , 

(28) 

where t 21 = r 2 — T\, r 21 = r 2 — rj , and zg = 2£m/j3 with 
I = 0, ±1, ■ • • . Using the Lehmann representation^ one 
can show that eq. (f26j) is connected with eq. (l28l) as 

p n ^(p,w) = -2Im£i,u(p, w + i0+), (29) 

where 0+ denotes an infinitesimal positive constant. On 
the other hand, £ can be expressed as eq. (|21a[) . whose 
poles are shared with those of /C in eq. (|21bj) through ■ 
Hence, the second mode also belongs to the two-particle 
Green's function K. Finally, x in eq. (|22b|) is not G 

itself due to the extra contribution -f( 2 ) + f (3) x (4) r (3) 
in the denominator. 

Thus, the two proofs by GSW may generally predict 
two distinct "massless" modes for the single-component 
BEC. As already mentioned, this fact has apparently 
been overlooked so far j 10 ' 11 it is certainly not incorpo- 
rated in the counting of Nambu-Goldstone bosons, which 
is based on the second proofs 

IV. PROPERTIES OF NAMBU-GOLDSTONE 
BOSONS 

Having clarified the differences between the two proofs 
of Goldstone's theorem by GSW^ we now study the 
properties of the corresponding two massless modes. 



X 

(a) 



(b) 



(O 



(d) 



(c) 



FIG. 2: Diagrammatic representation of the two-body inter- 
action in terms of the excitation field fa. A line and a filled 
circle denote fa and I^ - 1 , respectively. 
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(b) 



FIG. 3: Diagrammatic structures of {Hi nt ). A shaded circle 
denotes the proper contribution, which cannot be separated 
into two parts by cutting a single line. 



A. Improper structure of self-energy in BEC 

We start by clarifying a peculiar structure of the self- 
energy inherent in the BEC. Substituting eq. ((4]) into eq. 
(|2bj) . we can express Hi nt diagrammatically as shown in 
Fig. [5] in terms of the excitation field fa and the sym- 
metrized vertex of eq. (l24l) . Diagrams (b)-(e) are 
characteristic of the BEC. In particular, (b)-(d) give rise 
to the structures shown in Fig. [3] in the thermodynamic 
average (-Hint)- On the other hand, (H- mt ) can be ex- 
pressed analytically as eq. (jlOp . which was derived rigor- 
ously on the basis of eq. Comparing Fig. [3] with eq. 
(fTO]) . we are led to the conclusion that the self-energy it- 
self has the structure of Fig.[3][a), i.e., the structure called 
"improper" in the normal stated as it can be separated 
into two parts by cutting a single line. 

Improper contributions to 2 have been excluded in 
the classic studies of the BEC^r— i 12 i 14 ~ — where the self- 
energy appears to have been introduced by using the 
concept "proper" 3 ® without asking explicitly whether the 
relevant quantity is identical with the one in the Dyson- 
Beliaev equation. However, the above consideration has 
clarified that the distinction between "proper" and "im- 
proper" cannot be used to construct the self-energy de- 
fined by the Dyson-Beliaev equation (|Tb|) . 

The necessity of improper contributions to E was also 
pointed out recently on the basis of a self-consistent 
perturbation expansion simultaneously satisfying Gold- 
stone's theorem and conservation lawsi According to 
this expansion, the self-energies may be obtained from 
a single functional $ by eq. (|5aj) , and the functional <£> 
can be constructed order by order uniquely so as to sat- 
isfy eqs. © and (ITU|) simultaneously. Equations (|25a[) 
and (|25b[) present explicit expressions of the first- and 
second-order contributions, respectively. This functional 
4> also determines "irreducible" vertices in the Bethe- 
Salpeter equation (|20l) as eqs. (|18a[) - ()18dl) . In particular, 
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FIG. 4: Second-order diagrams for: (a) En and rj?; (b) 
£12 and . A line with an arrow (two arrows) repre- 
sents G (either F or F) in eq. (|23[) as in the theory of 
superconductivity.— The left number (the right number in 
round brackets) below each diagram indicates its relative 
weight in Eij (Iy ). Multiplying each number by 1/2 yields 
the absolute weight. 



eq. (|18dj) can be written in terms of $ as 

<5 2 $ 



ri?(i,2) = 4/3 



^(3,1)^3^(4,2) 
x(_l)j'-i+*+T 



%(3)*r(4) 



(30) 



This differentiation may also be carried out graphically 
by removing from the diagrams of $ a couple of lines with 
the second subscripts i and 3— j in all possible ways. Fig- 
ure [4] enumerates the second-order diagrams for Ey and 

{j = 1,2) obtained from eq. ([25b]) or Fig. □» by 
eqs. (l9al) and d30l) . respectively. Among them, the dia- 
grams in the bottom lines of Figs.HJa) andEJb) have the 
structure called "improper" in the normal state.™ Using 
the number below each diagram to indicate its relative 
weight, however, one can easily check that they cancel 
out exactly in the sum E + T^. Thus, they do not con- 
tribute to the denominator of eq. (|22b[) . The presence 
of improper diagrams in £ and their disappearance from 
E + r' 2 ) are expected to be a general feature of the BEC. 
We have also confirmed the statements in the third order 
of the self-consistent perturbation expansion using eqs. 
(Bl) and (B5) of Ref. G] and eqs. (JSaJ and (J30]) above; the 
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in the bottom lines of Figs. |U[a) and Hfb) in p space as 
(.7 = 1,2) 



£^ ip) (p) = 2g 2 n 2 [G 11 (p) - G 22 (p) - G 12 (p) + G 2 



(33) 

Now, if the retarded Green's function G R (p,uj) = 
G(p,uj + i0 + ) had a pole at w = E p , it would in 
turn bring about an infinite imaginary part at lo = E p 
in the retarded self-energies obtained from eq. (|33[) by 
zi — > lu + i0+. Thus, a sharp peak in the spectral func- 
tion 



A p (lu) = -2ImGii(p,w + iO + ) 



(34) 



FIG. 5: Various third-order diagrams for E and f (2) . (a) Di- 
agrams that contribute to both E and E + F 1 - 2 '. (b) Improper 
diagrams for E that disappear from E + r' 2 '. (c) Diagrams 
peculiar to f 1 ' 2 '. (d) Diagrams for I^ 2 - 1 with null total contri- 
bution. 



is impossible in the presence of improper contributions 
to £. In other words, no well-defined quasi-particle can 
exist in the single-particle channel. 

To see this more clearly, let us introduce an approxi- 
mation for G(p) where we incorporate only eq. (|33|) from 
the second-order along with the first-order self-energy 
studied in detail in Ref. l3ll This approximation 
may be quantitatively excellent for a dilute Bose gas; it 
will certainly enable us to qualitatively and explicitly sec 
how improper contributions affect the excitation spec- 
trum. Although eq. Q33p is apparently second-order in 

(2ip) 



i, i i • n • ttv Fn >ii i q, the self-consistency condition on Ei,,- will be seen 

results are summarized briefly m lug. [5] without arrows. V ■ j. -u • • 



B. Spectral function of a dilute Bose gas 

The improper structure of the self-energy £ inherent 
in the BEC has profound implications for the nature of 
its first mode determined as the poles of G. This mode, 
i.e., the Bogoliubov mode in the weak-coupling regimej 13 
has been regarded as a propagating mode with an infinite 
lifetime in the long- wavelength limit i 14 ' 16 ' 17 However, the 
improper structure gives rise to considerable broadening 
of the poles of G, changing the nature of the massless 
mode dramatically, as seen below. We call this mode, 
i.e., a mode with a substantial lifetime due to improper 
contributions to the self-energy, a "bubbling" mode based 
on the fact that it fluctuates temporally out of and back 
into the condensate reservoir. 

To convince ourselves of this statement, consider a ho- 
mogeneous system with the contact interaction 

V(ri - r 2 ) = gSfr - r 2 ), (31) 

and expand Green's function as 

(3(1,2) = J_^G(pK p -( ri - r2 )- z ^ Tl - T2 \ (32) 
v 



to yield a first-order contribution, changing the massless 
Bogoliubov spectrum dramatically. In contrast, other 
proper diagrams of the second order will only quantita- 
tively modify the first-order Bogoliubov spectrum to a 
tiny extent for the dilute Bose gas. Hence, they may be 
neglected for the present purpose. 

In general, eq. (I7b|) can be inverted for homogeneous 
systems to give 



Ep* - 



(35) 



where e p = p 2 /2m, E^ = £n(p), Ap = E i2 (p), Ep- = 
-E 22 (p) - £*^, Ap = -E 21 (p) = A*^, and Dp de- 
notes the determinant 



Dp 



•£ r 



Ep-- 



- f i) + ApAp. (36) 



The first-order self-energies E^ and A^ 1 ) have no p de- 
pendence for the contact potential of eq. (f3"T|) . Their ex- 
pressions are obtained by performing the differentiation 
of eq. (|9aj) with eq. (|25ap and carrying out the Fourier 
transform as 



£« = 2gn, A« 



p 



. (37) 



where V is the volume and p = (p, zi) with zg = 2£ni/ (3 
(£ = 0, ±1, • • • ). Noting that r<°) 2g and ^ ->■ ^ 
for the homogeneous system with no denoting the con- 
densate density, we can express the improper diagrams 



where 



n = no 



pv 



(38) 
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denotes the particle density. On the other hand, the ele- 
ments of eq. (|3"3)) satisfy 



■>(2ip). 



(2) 



(39) 



Adopting the approximation mentioned above, we can 
express the self-energies in eqs. (1331) and (133)) as 



Ejj = Ej, = + A< 2) , A p - = A = A« + A { ; 



Hence, eq. (fTaf reduces to 

M = S(i)-A«. 



(2) 

(40a) 



(40b) 



Using these results in eq. (|35]), we obtain a simplified 
expression for G(p) with Sjy — fi = Sp — fi = Ap — Ap. 
Let us substitute it into eq. (|33l) and use eq. (|59|) . We 
then find that eq. ([33]) forms a quadratic equation for 
Al 2) as 



with 



where 



(Af) 2 -( J D«/2 £p )Ag) + 2^=0 



4 1) -(^-^)(^+^ p b ; 



(£ P + 2AW) 



(41) 
(42) 

(43) 



is the Bogoliubov spectrum 7 - ' 17 ' 21 ' 31 with a linear disper- 
sion for |p| 0. By imposing A^ 2 ' — » as jz^j — ¥ oo, 17 



eq. (|4"Tj) is solved easily as 



><*) 



(i). 



1/2 



(44) 



Finally, the spectral function is obtained by eq. (jMJ), i.e., 



A p (ui) = -21m- 



£ P + A ? 



e p (s p + 2A P -) 



(45) 



which is a function of (ui, e p , A' 1 ), grto). 

It is worth pointing out that the spectral function (p4|) 
can be negative for condensed Bose systems. This may be 
realized by expressing it in the Lehmann representation 17 

A p (u) = 2tt^ ef( n -*"> [l - e-^'- £ ») 

x(ip^\c p \ip„>)((p„>\c p \ip^)S{uj -£ v , + £„), (46) 

where £„ and <p„ are an eigenvalue of the Hamil- 
tonian given by eq. (fTJ) and its cigenfunction, = 
— ln^ e~P £v is the thermodynamic potential, and 
c p and cp denote the annihilation and creation operators 
of the plane-wave state p, respectively. One may expect 
from eq. (|4"6"|) that E v i — £ v > holds for any finite ma- 
trix element (tf v i\c p \tf„) so that A p (uj) > 0. However, 
such an argument is no longer valid for condensed Bose 
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FIG. 6: Normalized spectral function j 4 p (oj)A' 1 ' for 
e p /A (1) = 0.01 as a function of w/E*. 

systems where cp is expressed as a linear combination 
of quasi-particle creation and annihilation operators.^ 7 - 
A classic counterexample is given by Bogoliubov theory, 



(2) 



whose spectral function is obtained by setting A^f 







in eq. (|4"5T) as 
A» = n[(a p + l)5(uj-E^-(a p -l)8(u 



-E*)}, (47) 



with a p = (e p + A^^/iJp 3 > 1. This function consists of 
a couple of sharp quasi-particle peaks at to = ±E p , and 
the one at u> = —E® is indeed negative. Note also that 
eq. (|47|) satisfies the general sum rule A p (u})dui = 2tt 
exactly, as it should. 

We have calculated eq. (|45|) by replacing gno — > A^ 
in the formula, which is justified for the dilute Bose gas 
near T = as seen from eq. (|37|). The resultant di- 
mensionless spectral function A p (w)A^ 1 ^ becomes a func- 
tion of only us/ A 1 ^ and e p / 'A^ 1 '. This approximation for 
A p (ui) will be excellent for weakly interacting Bose gases 
of gn/TQ <C 1 and A^ oc gn at low temperatures, where 
To = 3.31n 2 ' 3 /m denotes the transition temperature of 
the ideal Bose gas^i 

Figure [T5] shows A p (oj) as a function of ui for e p /A^ = 
0.01, which lies in the region of the linear Bogoliubov dis- 
persion in eq. (|43[) . Comparing Fig. [6] with eq. (j47j), we 
clearly observe that the couple of sharp quasi-particle 
peaks at ui/Ep = ±1 in Bogoliubov theory are sub- 
stantially blurred. Thus, the first mode is no longer a 
well-defined quasi-particle after incorporating the lowest- 
order improper contributions to the self-energy. Al- 
though the present consideration is based on the second- 
order analysis, the substantial lifetime is expected to be 
an essential consequence of the improper contribution, as 
may be justified by the general argument after eq. (|33[) 
or checked by investigating higher-order terms. Since the 
particle density n is given by eq. (|38l) . we may identify 
this "bubbling" mode as the temporal number-phase fluc- 
tuations of spontaneously broken U(l) symmetry. The 
absence of sharp poles in G(p) also has the effect of re- 
moving infrared divergences pointed out by Gavoret and 
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Nozieres. 4 Hence, the two-particle Green's function, for 
example, can be calculated safely without any artificial 
procedures by incorporating improper terms into the self- 
energy. 

Concerning the second mode, the improper contribu- 
tions disappear from the denominator of eq. (|22b[) . as 
already mentioned after eq. (1301) . Thus, wn l have 
a sharp quasi-particle peak in agreement with eq. (1271) . 
which may be observed as a propagating mode with den- 
sity fluctuations. We emphasize once again that this 
mode should be distinguished clearly from the Bogoli- 
ubov modeJ^ 



V. CONCLUDING REMARKS 

Perhaps it is not unreasonable to generally expect 
that single-particle and collective excitations of a many- 
particle system, determined as the poles of its single- 
particle and two-particle Green's functions, respectively, 
are different in character. Indeed, this is exactly the case 
for Fermi systems where the two excitations even obey 
different statistics of the Fermi and Bose types, respec- 
tively, and it is the single-particle excitation that is re- 
sponsible for most of their thermodynamic properties. 17 

The present paper has shown that the above state- 



ment will also hold true for the single-component BEC, 
contrary to the conventional viewpoint of attributing a 
unique phonon branch of collective density fluctuations 
as its excitation^— A characteristic of the system is that 
the two modes obey the same Bose statistics and may 
contribute equally to various phenomena. Thus, there 
remains a fundamental question: which mode dominates 
the thermodynamic, transport, and other properties of 
the single-component BEC? The answer may depend cru- 
cially on the interaction strength for the respective phe- 
nomena. In this context, it is worth pointing out that 
single-particle excitations are completely neglected in the 
variational wave function of Feynman for superfluid 4 Hej- 
which incorporates only the collective density fluctua- 
tions as may be justified in the strong-coupling region. 
However, single-particle excitations should also be rele- 
vant in the weak-coupling region as in the case of the 
ideal Bose gas. Further investigations are required to 
clarify this point. 
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